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Abstract In this paper, we review our main results involving the single particle momentum dis¬ 
tribution of bosonic trimer states in two and three dimensions. A summary table makes easier the 
comparison between the matrix elements and the different terms of the momentum distributions. We 
also show a practical method to continuously interpolate between different dimensions. 
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1 Introduction 

Many nature laws can be strongly affected if dimensionality is changed. As already pointed out by 
Landau in his classic book jj], any infinitesimal amount of attraction produce a bound state in 2D, 
while a finite amount of attraction is necessary to bind a 3D system. A remarkable phenomenon related 
to the dimensionality of the system arises in the study of three identical bosons, where the differences 
in the energy spectrum (and also other observables) are directly related to the number of dimensions 
that this system may access: in 2D there are only two three-body bound states linked to one two-body 
bound state in the limit where the range of the potential goes to zero On the other hand, in 3D, 
the number of three-body bound states may grow to infinity SI - this effect is now called by Efimov 
effect. 

The Efimov effect corresponds to an accumulation of the three-boson energy levels, toward zero 
energy, when the two-body scattering length tends to infinity. In this limit, where the two-body energy 
is zero, the energies of successive states are geometrically spaced obeying a universal ratio. These states 
were predicted and observed for three identical bosons in 3D systems |bj, but are absent in 2D even 
in the most favorable scenario of mass-imbalanced systems [Ji; 8], where a mass-dependent effective 
potential favors the binding of a light particle to a heavy dimer |9[. 

The appearance of Efimov states in 3D is very closely related to the possibility of collapse the 
three-body system. This collapse (Thomas collapse), firstly derived by Thomas in 1935, says that the 
three-body ground state energy may be made as deep as you want by decreasing the range of the 
potential (ro) - in the limit ro -*■ 0 the three-body binding energy tends to infinity. This divergence 
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demands the inclusion of a cutoff or, equivalently, a new physical scale independent of the two-body 
energy. In 2D, this collapse is absent in such a way the three-body observables are proportional to the 
two-body energy. For example, the three-body ground state energy is 16.52 E 2 and the energy of the 
first excited state is 1.267E 2 for three identical bosonsQ. 

There are many examples of observables in cold atomic gases that are affected by the dimensionality 
of the system. We would like to start mentioning the two- (C 2 ) and three-body (C 3 ) contact parameters. 
The connection between universal two-body correlations to many-body properties through the quantity 
C 2 was proposed by Tan in 0 (this quantity is often called Tan’s contact parameter). For example, 
the variation in the energy of a Fermi gas of momentum kp with the interaction strength (scattering 
length a) is directly proportional to this C 2l namely 

2n d\-l/(k F a)] = ° 2 ' (1) 


Furthermore, the virial theorem for this atomic gas also relates with C 2 through 


E- 2V = - 


C 2 

4-7T kpa 


( 2 ) 


These relations, in the way they are presented, were confirmed in experiments with two-component 
Fermi gases 0 , where each side of Eqs. © and © were measured independently and after compared 
to each other. A later experiment showed that similar relations also hold for bosons [l2j] . 

The quantities on the left-hand-side of Eqs. © and © are defined through the many-body proper¬ 
ties of the gas, while the contact parameter is defined in the few-body sector. A way to determine this 
parameter is to find the coefficient in the leading order of the asymptotic one-body large momentum 
density, n(q), of few-body systems, given by 


lirn n(q) -* + C 3 F(q) + ... . (3) 

q~>°o g 

The next order in this expansion defines the three-body contact parameter, C 3 , which may be im¬ 
portant only for bosonic systems, since the Pauli principle suppresses the short-range correlations for 
two-component Fermi gases. Notice that the momentum dependence of the leading order term in this 
expansion is the same for ID, 2D and 3D systems [l3|, but the function F(q) depends on the dimen¬ 
sionality of the system |14|. The two- and three-body contact parameters were determined for three 
identical bosons in 2D (l4l | and 3D and for mixed-species systems in 2D jl|| and 3D 0 - 

An interesting point about the two-body contact parameter is that, despite the considerable differ¬ 
ence between the binding energy of both states for three identical bosons in 2D (where the well-known 
limit cycle is not present), the ratio ^ is the same for the two states 0 - In general, for a mixed- 
species system - which have a richer energy spectrum !): Isf - the ratio is not the same for all 
states, but only in the special and experimentally accessible case of a three-body system composed for 
at least two identical non-interacting particles fl 6 j ]. 

Among the several differences involving the dimensionality of the system we would like to stress 
that the function F(q) in Eq. © has very distinct forms in each case of 2D or 3D. This function is 
directly related to the spectator functio ns f(q) (given in Eqs. © and ©), whose asymptotic form were 
discovered in the 60’s for 3D systems (ljjj and approximately 50 years later for 2D systems 0; fl 6 j . 
Thus, an interesting question is whether it is possible to interpolate between the 3D and 2D limits 
in a simple theoretical way and subsequently explore this in simulations using both more involved 
numerical methods and experimental setups since the development of the techniques for cooling and 
trap atoms allows the interpolation between different dimensions (13; (HI ■ 

We proposed a model that has the ability to interpolate geometrically between two and three spatial 
dimensions and thus study this crossover for both two- and three-body bound states of identical bosons. 
A “squeezed” dimension, whose size can be varied to interpolate the two limits, is employed with 
periodic boundary conditions (PBC). This model has the unique feature that it can be regularized 
analytically, which is a great advantage for its numerical implementation allowing to go smoothly 
between both limits. The theoretical elegance and tractability of calculations in the three-body system 
is itself a strong incentive for pursuing this geometry, but in spite of this elegance, a direct connection 
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between experiments and the parameter that dials between different dimensions with PBC in this 
model was not found yet. 

In the next sections we review our main results involving the momentum distributions in three and 
two dimensions. Both results are put together side-by-side in a table where the comparison becomes 
easier. In the last section we give an overview of our method that continuously interpolates between 3 
and 2D limits. 


2 Integral equation for bound states 

We investigate abc bound systems whose dynamics is restricted to either two (2D) or three spatial 
dimensions (3D). The masses are m a ,mb,m c and the pairwise interactions are described for attrac¬ 
tive zero-range potentials, being E a t, E ac , Ej, c the energy of each pair. The three-body wave function 
(q Q , p a^abc) has the same functional form in both 2D and 3D. For any s-wave bound state, the energy 
£3 is a solution of the free Schrodinger equation, except in the region where particles overlap. Using 
the Faddeev decomposition in momentum space, the bound state wave function in units of h = 1 is 
written as 

P' v mp+m-y ^ a |) f7 (|P a mp+m-i 

IF I 1 H 1 P* 

I h 2 mp-f 

where a, /3, 7 are cyclic permutations of a, b, c, q a is the a particle momenta with respect to the 
CM of the pair /3y, p a is the pair relative momenta, mp ia = m a (mp + m-y)/(m a + mp + m 7 ) and 
mp-y = (mp + m^Kjnp + m 7 ) are the reduced masses and f a ,p tl ( q) are the Faddeev components, or 
spectator functions. The three-body energy, £ 3 , and the spectator functions / Q ,/ 3, 7 are solution of a 
set of three coupled homogeneous integral equations, which in a compact form reads 

X oo 

k dk [K%p(q, k, £ 3 ) fp (k) + K^(q, k, £ 3 ) / 7 (k)] , (5) 

where the matrix elements of the two-body T-matrix, r^(q,E 3 ) and the kernels K®p(q, k, £ 3 ) and 
K® 7 (q,k,E 3 ) are given in table [ 1 ] 

An interesting difference between 2D and 3D three-body systems can be seen in Eq. © and table [TJ 
For each non-interacting /Fy pair, the respective spectator function f a {q) = 0. Choosing Ep 7 = 0 in 

the first line of table [Q gives different result for 2D and 3D systems. In 3D, [ 77 ^* (qr, £ 3 )] is finite 
and Eq. © is well-defined even if the three two-body subsystems interact with zero energy. On the 
other hand in 2D, [r^(q, £ 3 )] -*■ 00 , meaning that if at least two pairs have zero energy, Eq. © is 

not well-defined and three-body bound states do not exist. Therefore, non-interacting and zero-energy 
two-body systems lead to the same result in 2D, while they can give completely different results in 
3D [ 13 . 



fa(q a ) + fp(\ 

(q Q ,Pa Wabc) = ^ (q Q ,Pa) = - 


3 Momentum distribution 


The one-body density functions are observable quantities even in the limit of large momenta where 
the number of atoms is small, which has already been observed in experiments using time-of-flight and 
the mapping to momentum space (23j, Bragg spectroscopy 111 or momentum-resolved photo-emission 
spectroscopy fUJ. 

The one-body momentum density of the particle a is defined through the wave function ^(qa, Pa) 
from Eq. 0 as 

n(q a ) = f d°p a \<l , (q a ..pa)\ 2 ( 6 ) 


and the normalization is / d D q a n(q a ) = 1, where D = 2,3 for 2D or 3D systems, respectively. 

Inserting Eq. 0 in Eq. © and expanding it, the nine initial terms can be grouped into four 
components by using arguments of symmetry, each one with a distinctly different integrand structure. 
The one-body momentum density is expressed as a sum of this four terms, i.e., n(q a ) = n i(da)- 
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Table 1 Matrix elements of the two-body T-matrix, T E (q,Ez) and the kernels K®Aq, k, E 3 ) and (q, k, E 3 ) 
for both 2D and 3D systems, where p is the subtraction point (see, for example, [22111. 

2D 3D 


47rm^ 7 In ^ 

2 

9 Fz, 

2m fl 7 ,c E3 


>1 

( E 1 q2 1 fc2 ) ( k 1 1 

( B3+ 2m a7 + 2 m/37 J (m 7 ) 

1 

) 

\ 

( c. , I 2 . k 2 1 ( k q \ 

( E3+2 ™ a p + 2 ™-M) \ m p) 



A general system of three distinguishable particles, presents three distinct one-body momentum 
distributions, each one corresponding to a different particle. The four terms for particle a are expressed 
as 


ni(q a ) 
n 2 (q a ) 


n3 (q a ) 


n 4 (q a ) 


I f a (q a )\ 2 f 

/ 


d D k- 


\-E* 


k 2 


' l /3 7 ,o 


7 . 


d D k- 


\Mk)f 


( 


fa (da) J < 


2m Q 


k 2 

2771/3, 


Hi 

mry J 


I 


d D k- 


l/y(*)f 


\-E* 


2 m c 


k 2 


m p ) 


r d°k 

fe(k) 

m) 


( p, 1 I fc2 1 No) 2 

E\ , 9 2 . fc 2 k qA 2 


. y ^ 2m a -y 2mp*y m-y J 

' 2 m a p 1 2 mp 1 mp ) J 


/ 


d D k 




771 ft +771-y 


"•■I) 


( 


+ c.c. . 


-E 3 


p* 


2mp~f t cx 2mp*y 


( 7 ) 

( 8 ) 

( 9 ) 

( 10 ) 


Notice that the distributions for the other particles are obtained by cyclic permutations of (a,/3, 7 ) in 
these expressions. 

Although the equations for bound states (Eq. (0) and momentum distributions (Eqs. Q) to (fTOl) 1 
were derived for a general case of three distinguishable particles, we now specialize to experimentally 
relevant systems composed by two identical bosons a and a distinct particle b. The large momentum 
limit of Eqs. Q) to (flOl) were derived in detail for 2 D and 3D systems respectively in {T6j; LLZ1 and the 
final result for the density profile of particle b with respect to the pair aa is presented in table [ 2 j in 
units of m a = 1 . 

In the following we discuss two interesting properties in the momentum distribution of aab systems 
in both 2D and 3D. The geometric scaling of the Efimov states implies that observables may be 
described independently of the quantum state. The independence of the quantum states is not expected 
to be valid for 2D systems, since they do not present any geometric scaling. However, the leading order 
in the large momentum distribution was found to be independent of the state for three identical bosons 
[ 3 ; [H| and for aab systems, if the aa subsystem is not interacting |l 6 j. 

The effect of the two-body energy on the contact parameter is shown in Fig. |T] for the 40 K 40 K 6 Li 
system. This system has three excited states when E aa = E a t> and only two when E aa = 0. Notice 
that the large momentum limit of the momentum density goes to a constant in all cases. For E aa = 0 
both momentum distributions are equal in units of the three-body energy, i.e., n 0 (qb)/E 3 = ^(q^/E^, 
where the superscript denotes the quantum state. This case is rather special because the two identical 
particles have zero energy and cannot provide a scale such that the three-body structure is determined 
by the identical two-body interactions in the identical subsystems. In other words the large-momentum 
limit of the one-body density for particle a is determined by the properties of the ab subsystem. 
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Table 2 Asymptotic forms for both 2D and 3D systems where A = , tan 63 = for 0 < 63 < 7 t /2 and 

tan 04 = A(A + 2 ) for 0 < #4 < 7t/2, .T, Co and Cf, are normalization constants and n$(qb) is the second order 
term in the expansion of n 2 (qb) in Eq. ([H]>. 


ni(qb) 


2D 


167 r TT5 r: 


2 In ( 9 ;,) 


3D 



n 2 (qb) 


16tt .A 2 

If" (-A+ 1) 2 


/“dfcfc |/ a (fc )| 2 


32 7T A 2 

If (-A+1) 2 


/~dfcfc |/ a (fc )| 2 


n 3 (<76) 



n 4 (qb) 


87 rr 2 li^t) 

% 


n 5 (qb ) 


3277 ^4(-A-2) p 2 ln 3 (q 5 ) 
3 TXfTp" 


4-7T CgCfj 

'b coah (^f) 


■ cos I s In 


+ sin I s In 


sin ^ 


/ -4+1 
2.4 

/.a+t' 

2^4 


cosh [s (7 
sinh[s(f 


-®3)] 

•fa)]} 


877 2 |c a | 2 ^4 2 

s q 3 cosh(^) ^A(A+2) 


{v'A(^ + 2 )sinh[s(f - 6 ( 4 )] 
-j^ c ° s h[.s(f - 6 h)]} 


S7T 2 |e a | 2 -A 3 (.A+3) 

9? (.A+l)V-A(-A+2) 


This picture changes when E aa = E a b , as seen in Fig. [I] Now, in the large-moment urn limit, the 
coefficients of the one-body densities change with the excitation energy. The systematics is that the 
coefficients move towards the corresponding values for E aa = 0 as function of excitation energy. First 
the differences of the ratios with the two-body energies is understandable, since the interaction of the 
two identical particles now must affect the three-body structure at small distances, and hence at large 
momenta. However, as the three-body binding energy decreases, the size of the system increases and 
details of the short-distance structure becomes less important. 

The independence of the state in the one-body momentum distribution of 2D three-body system can 
happen or not, depending on the energy of the two-body subsystem. However, both mass-imbalanced 
and identical bosons systems present the same functional form to the leading order (LO) and next- 
to-leading order(NLO) in the large momentum expansion of the one-body density. The same does not 
happens in 3D, since the NLO contribution due to non-oscillatory terms vanishes when .4=0.20, 1.00 
and 1.57, i.e^ the sum of the components n\{qb),n 3 {qb) 1 n 4 : {qb) 1 nb{qb) from table [5] is null for these 
mass ratios (l7j |. 

We want to emphasize that both the results for 2D and 3D can be experimentally checked in the 
near future, since the calculations were made for alkali atoms. Besides the identical mass case, the 
other two mass ratios used in the 3D calculations can be obtained with mixtures of 133 Cs and 87 Rb 
{A = 1.565) or 6 Li and 39 K (A = 0.179). 


4 3D - 2D transition with PBC 

The physical and mathematical differences of three-body systems restricted to either 2D or 3D, pre¬ 
sented in the previous sections, are the motivation that lead us to pursuit a method where the dimen¬ 
sionality enters as a parameter allowing to continuously interpolate between the well-known extremes 
of 2D and 3D. We restrict analyses to three-identical bosons, which presents the Efimov effect in 3D, 
but only two three-body bound states in 2D. Furthermore, the dimensionality plays an important role 
in the momentum distribution already in this simplest case, as can be seen in table [ 2 ] for A = 1. 

Periodic boundary conditions (PBC) are assumed to be valid for the distance between the particles 
in the z direction. The relative momentum is given by p ± = (p x ,p y ) in the flat 2D surface and by 

\Z 2 nn n 
R ’ 


Pz = 


L 


ft — 0 , ± 1 , ± 2 , ■*•, 


( 11 ) 
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Fig. 1 The leading order term of the one-body momentum density divided by for each bound state labeled 
as n in a system composed of two identical (a = 40 K) particles and a distinct one (6 = 6 Li) as a function of the 
momentum q for both E aa = E a b and E aa = 0. 


in the transverse direction, with L = \Z2ttR being the size of the compact dimension corresponding 
to a radius R, which is the parameter that dials between two and three-dimensions. When R -*■ 0 it 
selects the 2D case and in the opposite limit, i.e., R-* o o, the 3D case is selected [26|. The momentum 
p and its corresponding phase factor dp are, with PBC, defined as 

n 2 1 

P = + ^2 and dp = Pl ' ( 12 ) 

We introduce the symbol f , which indicates an integration over the continuum momentum in the plane 
(p ± ) and a sum over the discrete perpendicular momentum (p z = ^). It reads 


p ^ l J 

J n=-oo J 



(13) 


Using dehnition m, the three-body free Hamiltonian becomes 


Hq (q- k ) = (q± + q -) 2 + ( k ± 




k z y 

„2 


(q± + q z ) ■ ( k i 

„2 ^ ^ 


k z) , 


" 1 + qi ' k 1 + & + lv + lv 


(14) 


and considering Eqs. (HI to m, the integral equation for the bound state ([5]) for a compact dimension 
with PBC is found to be 


f(<l±,n) 




/(ki,m) 
—E 3 + Hq (q, k) 


/ (k±, m) \ 
P 2 + H p 0 ( q,k)j 


(15) 


with r p (E) given in Eq. (fTCl) and Hq (q, k) in Eq. d. The subtraction is kept even after the dis¬ 
cretization because the Thomas collapse is always present for any finite compact radius, no matter 
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how small it is. It is worthwhile to remind that, for R -*■ oo, Eq. USD returns precisely the equation for 
the spectator function in 3D 

The two-body scattering amplitude is 


r P {EY l 


2t t sinh^Tr^yjEf) 

R sinh 


(16) 


which recovers the matrix elements of 3D and 2D systems in the limits R -> oo and R -> 0, respectively. 
The first case is straightforward an reads 


7-3“ b(E) = = -2tt 2 (v^i- V\e 7 \) . 


(17) 


Going to the 2D limit, it is important to notice that a quasi-2D system is in practice a 3D system. 
Then, the units of t^(E) and r p 1 {E) are exactly the same, but are different from t^^E). Taking 
into account the correct units, the 2D limit of Eq. ED reads 


r 2 d(E) = lim R r p l (E) = -27rln 
R—> 0 



(18) 


Expressions in Eqs. (ED and (1181) are respectively identical to the expressions presented in table [T| for 
3D and 2D two-body T-matrix when m a = nib = m c - 

Introducing dimensionless variables, £3 = F 3 // 1 2 , £2 = F 2 / y 2 , r = R n, y ± = q±/^JI and x L = k ± /^/Jl 
and integrating over the angular dependence, since the focus is on states with zero angular momentum, 
the integral equation ED is written as 




7rln 


sinh (nryj'±(yl + (£) -e 3 j 

sinh (-7T r>/ej) 


-1 


00 c 

Z / dx ± x ± --- 

0 \\f(-e 3 + yl + xl + £ + 1 £ + ^) -ylxl 


^J(l + yl +x l + n? + !n? + ^) 2 - x l yl 


■ (19) 


The dimensional crossover transition is explored through the numerical solution of Eq. (flTH) . In 
Fig. [2] the ratios £ 3/^2 are showed as function of the compact dimension radius r, for the ground, first, 
and second excited states. Notice that the last state goes into the continuum before the 2D limit is 
reached. 

The computations were performed for two fixed two-body energies £2 = 10 ~ 6 (empty circles/dashed 
lines) and 10 ' (full circles/solid lines). Note that the Efimov ratio between two consecutive three- 
body states, ~ 515, is not completely reproduced for a finite a. The points at which the energies are 
calculated are showed explicitly, while the curves are guides to the eye and for r = 1000 the energies 
are obtained from the pure 3D equation. 

An interesting dimensional crossover result is seen in Fig. [2} where only one sharp transition is 
present for the ground state while there are two for the first excited state. This behavior can be 
understood by considering the size of the trirner given roughly by r ~ 1 IsfEj,- For £2 = 10 7 , the ground 
state plateau for £ 3/62 = 93330 is placed at f = 10.35 and first excited state plateau for = 211.79 
at f = 217.29. These f values give approximately the region of the jumps signaling that the 3D limit, 
represented by the plateau, is reached once the trirner size matches the size of the squeezed dimension, 
r. The same analysis can be made for £2 = 10 ~ 6 with f = 10.27 and r = 188.98, respectively, for the 
ground and first excited state. Varying r from large to small values, the 3D->2D transition occurs 
for r ~ 10 , where it is possible to notice the disappearance of the higher excited states in order to 




















r 

Fig. 2 63/62 as a function of r, for 62 = 1CT 7 (full circles) and 1CT 6 (empty circles). The solid and dashed lines 
are guides to the eye. As the 2D limit (r -> 0) is approached, higher excited states disappear and only the 
ground and first excited states remain. 


reproduce the well known 2D results with two trimer bound state energies proportional to £2 with the 
ratios £ 3/62 = 16.52 and e 3 /e 2 = 1.27 @]. 

From the experimental point of view it may be difficult to keep the dimer energy constant. However, 
the transition observed in Fig. [2] will not disappear due to a variation of e 2 with r. The increase of 
the dimer energy will merely move the beginning of the jumps towards smaller r. The optimal way 
to probe these jumps is to start from a two-body energy in the unitary limit (a -*■ 00 ) where the 2D 
plateaus are fixed. Larger dimer energies will cause the 3D plateau to move to lower £ 3/62 ratio and 
push the beginning of the transition to smaller r, thus making the transition region broader. Another 
interesting study about the dimensional crossover, where three identical bosons are confined by a 
harmonic potential along one direction is found in |27l |. 


5 Conclusion 

In this paper we summarized our main results involving the single particle momentum distributions 
in two and three dimensions. The summary tables, comparing the matrix elements and the functional 
form for the different terms of the momentum distribution makes the comparison between the 2D 
and 3D regimes much easier. In the last section we presented in a very schematical form a method to 
continuously interpolate between different dimensions. The same technique used here to go from 3D to 
2D, may be used to go from 2D to ID systems. All discussions in this text followed closely the papers 
jl6c I17t 1261 ]. An interesting direction for future investigation, which connects all the results discussed 
in this work, is to understand whether the contact parameter tell us how much of the wave function is 
in each dimension in the transition region presented in Fig. 10. 
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